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Recently Bar-Ziv and Moses discovered a dynamical shape transformation induced in 
cylindrical lipid bilayer vesicles by the action of laser tweezers. We develop a hydrodynamic 
theory of fluid bilayers in interaction with the surrounding water and argue that the effect 
of the laser is to induce a sudden tension in the membrane. We refine our previous analysis 
to account for the fact that the shape transformation is not uniform but propagates out- 
ward from the laser trap. Applying the marginal stability criterion to this situation gives 
us an improved prediction for the selected initial wavelength and a new prediction for the 
^ . propagation velocity, both in rough agreement with the experimental values. For example, 
a tubule of initial radius 0.7/um has a predicted initial sinusoidal perturbation in its diam- 
eter with wavelength 5.5/im, as observed. The perturbation propagates as a front with the 
qualitatively correct front velocity a bit less than 100/um/sec. In particular we show why 
this velocity is initially constant, as observed, and so much smaller than the natural scale 
set by the tension. We also predict that the front velocity should increase linearly with 
laser power. Finally we introduce an approximate hydrodynamic model applicable to the 
fully nonlinear regime. This model exhibits propagating fronts as well as fully-developed 
"pearled" vesicles similar to those seen in the experiments. 
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1. Introduction and Summary 

The study of artificial biomembranes has opened a window into the machinery of 
real cells by giving us physical systems which are simple enough to describe from first 
principles, yet complicated enough to display lifelike behavior. The study of the equilibrium 
configurations of closed bags of lipid bilayers ( "vesicles" ) is by now well advanced |l| . 
Most interesting biophysical phenomena, however, are not in equilibrium, and the study 
of the dynamics of vesicle shapes is somewhat less developed. For example, budding and 
vesiculation || and instabilities crucial for understanding adhesion ||, are all inherently 
dynamical processes. 

To gain a physical understanding of non-equilibrium membrane dynamics we must be- 
gin with experiments in which some physical intervention brings about a dynamical shape 
transformation.0 Several such techniques are by now well developed. Budding and other 
shape transformations can be induced by adjusting temperature in closed vesicles [f|[||. 
Strings of beads can be induced by relieving tension using micropipettes 0, by sudden 
hydration of dry lipid || , and by the incorporation of a small object which is subsequently 
pulled away with laser tweezers [0 . One can even create tension inside a vesicle by poly- 
merizing a stiff rod inside it [!(]] , whereupon a cylindrical structure modulated with pearls 
can emerge ||1 1|| . 

Recently Moses and Bar-Ziv have introduced a new class of experiments in which 
laser tweezers act directly on single lipid bilayers JT^]. It is already perhaps surprising 
that there should be any action at all, since the optical absorption of a single bilayer is so 
small. As we will recall in Sect. 3 below, we argued in |13[ that the effect of the tweezers 
is to induce a tension E in the membrane proportional to the laser power. Our argument 
relied only on basic electrodynamics, and while it does not give the precise value of the 
constant of proportionality, it does imply that the effect of the laser is extremely simple; 
unlike other more mechanical probes, the possibility of complicated parasitic effects seems 
minimal. More importantly, laser tweezers provide for the first time the fast, delicate, and 
highly localized probe needed to understand membrane dynamics in detail. For example, 
while beaded tubes have been seen in some of the other techniques mentioned above (see 
Sect. 2), a unique feature of the tweezer experiment is the controlled ability to excite a 
"small-amplitude" (quasilinear) peristaltic modulation. 



1 Experiments in which a biological intervention is introduced (for example, mutating a 
single protein) are comparatively abundant 0. 
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We recall the basic results of [|12] in Sect. 2 below. The most intriguing phenomenon for 
us was the fact that when laser tweezers are applied to previously stable long cylindrical 
vesicles, they excite an instability ending with a "string of pearls" state. The initial 
wavelength of the instability is A = 2nR Q /k, where -R is the initial cylinder radius and the 
dimensionless wavenumber is typically k = 0.8. Remarkably, a purely local excitation due 
to the laser trap creates (in a certain regime) a uniformly modulated state, which invades 
the initial cylindrical region at a roughly constant velocity of about Vj- ~30 / um/sec. This 
is not the sort of behavior we usually associate with pulling a stretchy object in a viscous 
fluid, and so we have a very interesting dynamical system. 

The pearling phenomenon bears a superficial resemblance to the instability of a cylin- 
drical column of water in air, studied in the classical works of Plateau and Rayleigh 
||14j| [I5| PB| . In this situation surface tension destabilizes the cylinder, since the same 
volume per unit length can be contained with less surface string of spherical 

droplets. As in other pattern-forming systems, a competition ensues: the lowest surface- 
to- volume ratio comes from a small number of large droplets, but the formation of such a 
state is kinetically suppressed because it requires the motion of water (a conserved quan- 
tity) over long distances. For a macroscopic system the water may be treated as inviscid 
(high Reynolds number); Rayleigh calculated that the fastest-growing instability set in at 
wavenumber k = 0.70 [|15"f . At micron scales, however, we are in a regime of low Reynolds 
number; here Rayleigh found the fastest-growing instability to be at k = ||16|1 , totally 



different from what is seen in [12 



Of course the situation studied in ref. |T^] is not a thread of liquid surrounded by air. 
Long ago Tomotika considered a two-fluid model, in which a column of one viscous fluid is 



initially immersed in another with a certain positive surface tension |17j . One could imagine 
that these two fluids correspond to the interior and exterior water, so that both viscosities 
are equal. Specializing Tomotika's result to this case, one finds the fastest-growing mode 
to be at k = 0.56. Unfortunately this simple picture is not obviously adequate for the 
pearling problem, nor indeed is its predicted value for & max ever observed; experimentally 
the initial wavenumber is always a bit larger, and then increases still further with time. 
What is missing? 

For one thing, the two-fluid model neglects the presence of a material object between 
the inner and outer fluids, namely the membrane, except insofar as the latter somehow 
communicates tension induced at the trap to the entire surface. In particular Tomotika 
imposed continuity of the shear stress T pz across the boundary, while a material object 



located there will in general be capable of exerting tangential forces and hence causing a 
discontinuity in T . 

A more realistic treatment of membrane dynamics requires a three-fluid model, in 
which a two-dimensional fluid (the bilayer) separates two three-dimensional fluids (the 
water). The "tension" on the membrane can then more precisely be regarded as minus 
the 2d pressure IT of the intermediate fluid, which in turn emerges from the appropriate 
hydrodynamic equations. The key early works on membrane dynamics used models of 
this type to study equilibrium fluctuations jlBJ . In the pearling experiment the laser takes 
the system far from equilibrium by suddenly imposing a boundary condition on IT. In 
other words, the laser trap is regarded as a lipid reservoir located at z = in cylindrical 
coordinates, whose chemical potential suddenly jumps from zero to a negative value as the 
laser is switched on. 

The idea that the laser effectively induces a tension E, and that a Ray leigh- type 



instability ensues, was first proposed by Moses and Bar-Ziv |]I2|| . We substantiated this 
picture of the laser action in [|13[j , showed how to estimate E from the laser power, and 
found the fastest-growing mode /c max = 0.68, at the low end of the observed values.! 

The most serious lacunae in the analysis of ref. [ITBI were that (a) we neglected the 
propagating character of the pearling instability, assuming instead that the tension was 
everywhere a constant, and (b) we neglected the fact that the lipid fluid is conserved. 
These points are related, of course. For a shape transformation to propagate outward 
from the trap, leaving behind a stationary shape, lipid must constantly be transported 
away from the moving front and into the trap. Assumption (a) gave us no possibility at 
all of predicting the front velocity, a readily accessible experimental quantity. 

In this paper we will present a more detailed analysis addressing these and many 
other points. After reviewing the experimental facts in Sect. 2, we will build up a physical 
picture of the pearling system in Sect. 3. In particular we will show that assumption (b) 
above is a good approximation during the time regime of the experiment, so that pearling 
is essentially the invasion of a saturating front into a uniform linearly-unstable state. To 
get to the point as quickly as possible we will construct in Sect. 4 a very simple form of the 
three-fluid model which captures most of the essential physics without lengthy formulas. 
We then analyze the front propagation in Sect 5 using the marginal stability criterion (or 



2 Actually we quoted /c max = 0.65 due to an error which we will correct in eqn. 4.2 
below. 
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"MSC") ||19|| [20|j and show how it gives qualitatively correct results for and k when 
applied to our simple model. As we will recall, the MSC allows one to compute the front 
velocity using only the linearized dynamics. 

To get quantitatively accurate predictions we will construct a more realistic model in 
Sect. 6. The reader may wish to skip this rather technical section. In it we first solve the 
hydrodynamics exactly, then resolve the bilayer into its two monolayers, to get a four-fluid 
model. The motivation for this is that dynamical friction between the two leaves of the 
bilayer has been shown to be quite significant for a related class of membrane problems 
[21}] ||22|| p3| , though as we will see the actual effect on our answers will be slight. Again 



using the MSC, we first find that the initial wavenumber is k = 0.80, which agrees with 



experiment somewhat better than the result of |T3|]. In particular the initial wavenumber 
is insensitive to changes in the laser power or tubule radius, as observed. Secondly we get 
a front velocity v j- equal to O.O6E/77 where 77 is the viscosity of water, another prediction in 
rough agreement with experiment. While this prediction isn't precise due to our imprecise 
knowledge of E, it does explain why the magnitude of v j is observed to be so much smaller 
than the natural velocity scale S/77. Moreover, the linear dependence of the initial front 
velocity on the laser power should be verifiable in the future. 

Since the MSC uses only the linearized equations it must assume, rather than prov- 
ing, the existence of a uniformly-propagating front. One could take the phenomenological 
attitude that such behavior is observed experimentally, but this and other qualitative facts 
should emerge from the solution of the full nonlinear equations. In Sect. 7 we introduce a 
third model, with most of the simplifications of the first, simple model (Sect. 4) but retain- 
ing the nonlinear structure of the elasticity. We show the results of numerical studies that 
indicate that it indeed supports propagating front solutions, and that it produces station- 
ary pearled structures like those seen in experiment. More details will appear elsewhere. 

Appendix A contains a glossary of symbols used for physical constants, and their 
values. 

Granek and Olami were the first to attempt a realistic treatment of pearling without 
the simplifying assumptions (a), (b) above |2|]. While our physical picture and conclusions 
are different from theirs, we are indebted to them for emphasizing these issues. We are 
also grateful to Moses and Bar-Ziv, who first suggested to us that marginal stability might 
be applicable to this problem. Finally, others have suggested that the laser could have 
other effects besides inducing tension, for example effectively inducing a bilayer asymmetry 
(spontaneous curvature) f24]l |]25|1 . We don't see how this could happen, and in any case we 
will see that no such effect is needed in order to explain the observed phenomena. 
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2. Relevant Experimental Facts 



Briefly the observed phenomena of interest to us are as follows [12| . Initial preparation 
of the system yields stable, nearly straight cylinders up to hundreds of microns long, 
anchored at both ends by large globules of lipid. Each tube consists of a single bilayer 
of DMPC or DGDG, with water on the inside and outside. The tubules are polydisperse, 
with initial radii R Q between 0.3-5/um. The high temperature used (45°C) implies that 
the membrane is in its pure fluid state. Initially the system is somewhat flaccid, as seen 
from visible thermal undulations and the fact that the tubes are not quite straight. 

Application of a laser spot localized to ~ 0.3/im produces a dramatic transformation 
to a stationary "peristaltic" figure, i.e. a cylindrical shape with radius r(z) at first varying 
roughly sinusoidally with distance z from the trap. There is no evidence for any initial time 
delay between laser illumination and the onset of the instability. Greater laser power is re- 
quired for larger-radius tubules, but nothing seems to depend on the length of the tube, so 
long as the trap is initially many radii from the ends. The shape transformation propagates 
outward from the laser trap, with a well-defined velocity Vf which varies between exper- 
imental trials but is typically about 30^im/sec and roughly constant for at least a couple 
dozen wavelengths. Remarkably, after a very short illumination the shape transformation 
continues to propagate after the laser is shut off, leading to a uniform, small-amplitude 
peristaltic shape. Longer excitation leads to the eponymous pearled state. Sometimes 
tubes intersect each other; in this case the shape transformation can actually cross from 
one tube to the other upon reaching the intersection ||26[| . 

Once formed, the peristaltic shape has a well-defined initial wavelength A which is 
uniform over many microns. The dimensionless initial wavenumber k = 2nR Q / A is always 
found in the range 0.64-1.0, and is typically about 0.8 (ref. |T!J, fig. 1), whatever the initial 
radius R . After prolonged tweezing some buildup of lipid becomes visible at the point 
of application of the laser. As the modulation grows more pronounced, k grows from k 
to become slightly greater than 1 and deviations from a simple sinusoidal profile become 
pronounced. The modulated state is tense: visible thermal fluctuations are suppressed and 
the tube draws itself straighter than initially. 

We pause to contrast the above phenomena with other related results. Deuling and 
Helfrich studied aged red blood cells, which had different interior and exterior fluids; in 
this situation it is reasonable to invoke a spontaneous curvature c expressing the chemical 
asymmetry of the bilayer environment. In our case there is no such asymmetry; as we will 
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recall in Sect. 3.1 below, the fact that each leaf of the bilayer is initially in equilibrium 
with a common lipid reservoir then implies that c = |3yJ . 

Evans and Rawicz, and recently Pouligny, have rapidly formed thin tethers coming 
from large vesicles, respectively by adjusting the internal pressure or by pulling out a 
small inclusion in the membrane || . Such tethers sometimes contain pearls. Also lower- 
temperature experiments, where in-plane order is important, can give pearled shapes PS]. 
All of these experiments have in common that the formation of pearls is coterminous 
with the formation of the tethers themselves, and small-amplitude sinusoidal modulations 
are not seen. Finally, experiments in mixed surfactant systems have shown equilibrium 
peristaltic shapes [^9[ . This is presumably due to shape-composition coupling, which is 
known to lead to an instability towards an asymmetric bilayer | 3"P|1 |3~IP and thence to the 
mechanism of 



3. Physical Picture 

In Rayleigh's problem, the tension was a material parameter characterizing the air- 
water interface, and hence trivially constant. In our case we argued in the Introduction 
that tension is instead a dynamical variable, and so its spread must be self-consistently 
determined along with the change of shape. Our main objective in this section is to justify 
the assumption [[T2J ]T3[] that tension initially spreads so rapidly that it effectively becomes 
constant and uniform as soon as the laser is turned oni. Thus the spread of tension is 
not what limits the speed of initial propagation of the shape pulse. Our arguments in 
this section will all be rather crude. We can only emphasize that the field of membrane 
dynamics is not yet fully developed, and buttress our arguments whenever possible with 
the observed phenomenology. 

In this section we will for simplicity assume axial symmetry everywhere. We will 
return to this point in Sect. 6.3. Thus all scalars are functions of radius p, distance z from 
the trap, and time t\ all velocity vectors have vanishing azimuthal component. 



As mentioned earlier, Granek and Olami first studied this issue 24 



3.1. Initial equilibrium 



We begin with a discussion of the initial equilibrium, before turning on the laser. 

Initially some preparation protocol has created long cylindrical vesicles constrained 
to stretch between the terminal blobs and in thermal equilibrium with them. We will 
speculate as little as possible about the nature of the blobs. Since they are far from the 
illumination spot, their only role is to determine the nature of the initial equilibrium, 
not the subsequent fast dynamics. We will simply assume that the blobs furnish lipid 
reservoirs and that both leaves of the bilayer membrane are in equilibrium with these 
common reservoirs, with a chemical potential for exchange close to zero. 

With this assumption we find that initially the membrane has tension close to zero, 
consistent with the initially observed thermal motion |12| . We also get that the spontaneous 
curvature c = 0, since in a chemically symmetric bilayer c can only arise as the difference 
in chemical potentials between the two constituent monolayers^ ||27|| . c = is consistent 
with the observed absence of any preferred initial tubule radius [THl . Nonzero spontaneous 



curvature can lead to equilibrium unduloid shapes fl32fl|[33| , but as we mentioned earlier 
these are not observed in the experiments in question. Moreover, these "Delaunay surfaces" 
begin with an initial wavenumber k = 1, larger than what is observed (see Appendix B). 

The physical boundary condition imposed by the terminal blobs will not matter in 
our analysis, but for concreteness suppose that they seal off the tubule, imposing fixed 
volume. This constraint then creates a physical pressure difference, with a very small 
negative constant pressure inside the tube p = —k/2R 3 (and p = outside), where k is 
the usual bending modulus. p balances the tendency of a cylinder to increase in diameter 
to reduce its bending energy. One can either solve the volume constraint explicitly and 
substitute into the bending energy (as in £l3[]), or else regard p as a Lagrange multiplier 
(as in ||33|| ) to see that then the cylindrical shape is stable to all small perturbations. This 
is the initial equilibrium state. 



4 See Appendix E. The situation is quite different in closed systems, such as vesicles of 
spherical topology. Here indeed the initial preparation leads to a new parameter describing 
vesicles, the fixed area difference between the two monolayer leaves || . 



3.2. Laser action 



Next we recall our model of the laser action from ref. |fL3l. When the laser comes 



close to the membrane, nothing happens: local heating is not important. When the laser 
spot touches the membrane, it pulls material in by the dielectric effect. While it is hard to 
calculate the exact tension so induced, we may easily estimate it as followsi: the applied 
laser power of ~ 50mW, spread over a spot of diameter 0.3/um, corresponds to an energy 
density in vacuum £ of 2.3 • 10 erg cm -3 . Taking the dielectric contrast between water and 
lipid at this frequency to be of order 5e = 0.23 [S3], we see that when a lipid molecule falls 
into the trap displacing water we gain an energy |SJJ ~ £5e ■ a D, where a is the area of 
the molecule's head and D is its total length. Taking 2D ~ 40A, we get that each unit of 
bilayer area sucked into the trap yields an energy gain of £ ~ 2 ■ 10 -3 ergcm -2 . 

Actually this figure is surely an overestimate. If the action of the laser is really to 
pack more lipid into the trap, there must be an offsetting cost per unit area to fold it up or 
otherwise put it into a more compact configuration than a single nearly-flat bilayer. (For 
more details see jSJ|.) Still, we see that the trap generates a tension well in excess of the 
critical value [|I2] S crit ~ ~ 1.2 • 10 _4 ergcm -2 needed to trigger shape transformations 
(see Sect. 6.3 below), where k ~ 0.6 • 10 _12 erg is the bending stiffness of DMPC bilayers 
and we took R = 0.7/um for illustration. Let a = T,R 2 /k denote the dimensionless 
tension; thus 

a = ER q 2 /k ~ 2 • 10 1 , (3.1) 

and larger for larger tubules. Such large tensions are qualitatively corroborated by other 
experiments involving laser-induced expulsion of vesicles [36[. 



One could imagine that the details of how the bilayer gets packed into the optical 
trap could effectively create different boundary conditions for the chemical potential of 
each layer and hence induce a sort of symmetry breaking, a spontaneous curvature for the 
layer. For example, tiny vesicles could get pinched off preferentially from the inner layer. 
We will resist such speculations, as we have no theoretical or experimental evidence in 
favor of such a mechanism, nor will such a spontaneous curvature be needed to understand 
the phenomena. We will also assume that the laser creates no major disruption to the 
integrity of the bilayer, i.e. the membrane remains impermeable on the time scale (tenths 



We thank R. Bruinsma for suggesting this estimate. 
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of a second) in question!, and transfer of lipid from one leaf to the other ( "flipflop" ) also 



remains too slow to be of interest, as in unstressed membranes |37[ . 

Further support for our model of the origin of pearling by laser-induced tension comes 
from the observation that mechanical tension applied suddenly to a cylindrical vesicle (by 
dragging an end attached to a movable pipette) creates a pearling instability similar to 
the one seen with the laser |26|| . Whichever excitation one uses, when it stops the tension 
reverts slowly to zero, thermal fluctuations resume, and the tubule relaxes back to its 
initially stable cylindrical shape |12 . 



3.3. Initial tension propagation 

As we mentioned in the Introduction, the induced tension E is properly to be regarded 
as a boundary condition on the pressure IT of the lipid fluid@: 

U(z = 0,t > 0) = -E . 

The sudden introduction of a gradient of IT causes the lipid to move toward the laser spot; 
this motion in turn stretches the membrane, causing the tension to spread outward. We 
will first argue that this initial spreading of tension is very fast, and essentially complete 
long before shape changes have had a chance to begin. Hence in this subsection only let 
us at first neglect shape changes altogether and ask what happens in a linearized analysis 
when we begin to pull on a cylinder of 2d fluid, then later check the self-consistency of this 
picture. Since the shape is fixed and the water is incompressible, there can be no net flow 
of water across a cross-section of the tube. 

As the bilayer begins to move toward the trap, it entrains the surrounding water. 
Focusing our attention on the interior volume, this entrainment sets up a velocity profile 
v z (p, z,t). For this estimate we will use Poiseuille (or "lubrication") approximation, where 
all gradients in z are assumed much smaller than those in p. (We will also temporarily 
neglect the exterior fluid.) The incompressible conservation law V • v = then implies 



6 This amounts to assuming that the volume is effectively clamped by dissolved macro- 
molecules. It is reasonable since large laser-induced tension can pressurize vesicles with 
no significant loss of interior volume, even long after the laser is shut off [[J6|. Indeed far 
greater tensions than those considered here can be applied mechanically with no observed 
loss of interior volume nor other breakdown in the bilayer model 0. 

7 The pressure of a 2d fluid has dimensions of tension. 
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that the radial velocity is small, v <C v z , and by the equations of motion the pressure p 
is constant across the cross-section of the tube. Solving the remaining equation of motion 
rf\/ 2 v z = V z p in this approximation then gives us the usual parabolic velocity profile of 
Poiseuille flow, 

v z (j>,z,t) = f(z,t)p 2 -g(z,t) , (3.2) 

where / is related to the pressure gradient and g = \Rq f because there can be no net 
flow of water down the tube. 

The membrane itself is a compressible 2d fluid. Its velocity v must match that of the 
water: 

v = v z (R ) = \R 2 f 



For its equation of motion, we can neglect the 2d viscosity of the lipid on these scales |22| 
and write the force-balance equation as 

V Z II(M) = T zp (R Q ,z,t) = -vd p v z \ Ro = -2riR f(z,t) , 

where is the 3d stress tensor of water. Gradients in the velocity must in turn affect the 
density of lipid, by another conservation law. We write the lipid density as <$> = O (1 + x) 
where O is the equilibrium density and x is the areal strain; then 

— [dz 2ttR <P (1 + X )] = -dz- V z [2ttR v ■ O (1 + X )} ■ (3-3) 

To close the equations we need the constitutive relation 

n = ^x , 

where K is the 2d bulk modulus of the lipid layer. 
Combining we find 

dU KR n 2 , A . 

whose solution H(z, t) indeed spreads rapidly until it is essentially equal to the boundary 
condition — E throughout the observed tens of microns. To estimate how rapidly, we need 
the value of K. We argue in Appendix C that K is effectively much smaller than its "bare" 
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value! K Q ~ 1.4 • 10 2 ergcm 2 ; instead we will argue for K eS ~ 10 1 ergcm 2 . Even so, 
K eS ^> E. Combined with (|3_7ji]) we § e ^ the hierarchy of scales 

kR ~ 2 < E < K eff , (3.5) 

which we will use repeatedly. 

The large modulus means that the tension profile will rapidly outrun any front of 
shape transformation traveling at v* < E/77 (in fact we will find in Sect. 6 that v j is much 
smaller than this). For technical reasons the observed region must contain the laser spot, 
so propagation is not observed for values of L larger than about 50/um. Scaling ([J74]) we see 
that at a distance L from the trap the tension approaches its saturation value in a time of 
order t L ~ L 2 r)/K eE R Q . Taking L = 50/um gives t L ~ 4 • 10~ 2 sec, about one video frame, 
and much faster than the time L/vi ~ 2 sec needed for the observed front propagation. 

We conclude that after a very short time the tweezers create a uniform, tense, stretched 
state of the lipid in the region of interest. If the tension E is great enough, this state will 
be everywhere linearly unstable to shape perturbations, just as in the Rayleigh instability. 
Of course it remains to show that a uniform stretched state will indeed change shape via 
a uniformly propagating front as assumed above. 

3.4- Front propagation 

We have just argued that initially the tension spreads rapidly, before the shape has 
had a chance to change. To see what happens to tension as the instability progresses, we 
will continue to assume that this proceeds by the propagation of a front, leaving behind a 
stationary modulated state. Later, this is shown to be a self-consistent picture. 

The energy which drives our instability comes from delivering lipid to the trap. But 
where does this lipid come from? Certainly very little comes from stretching the bilayer, 
since we just asserted that its bulk modulus is much greater than the applied tension. 
Instead it comes from a shape transformation, in which the initial cylindrical shape gets 
replaced by something with less area (and hence fewer lipid molecules) per length. It may 
at first be hard to see how such a transformation could propagate at a constant velocity. 
After all, lipid must constantly be transported to the trap, entraining water as in the 



8 Were we to use the bare value K , we would have to include inertial effects in the 



above derivation, leading to a different dispersion relation from (37 4) but with the 



same qualitative conclusion: tension spreads more quickly than shape change. 
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previous subsection and incurring viscous dissipation. Much as in Poiseuille flow through 
a pipe, shouldn't we expect a decrease of pressure |II| as the front moves away from the 
trap? 

To address the question, consider first a case in which the cylinder gets converted at 
a moving front to a modulated shape, with a change of its area per unit length from 2nR 
to 2ttR (1 — a). The front is located at z = L, so to advance it at a velocity = dL/dt 
lipid must flow toward the trap at velocity v = aVf. 

The laser does work on the membrane at a rate 2'kRqOlv^ ■ E. This energy is lost due 
to viscous dissipation both at the front, and everywhere from z = to L. Since the front is 
only a few times R in length, we can estimate the former loss as r]Vf 2 7:R Q . The latter loss 
will be mainly due to entrainment of water. By the remarks in the previous subsection, 
the entrainment of water creates a shear d v z ~ avj/R in a volume tcR 2 L, for a loss 
rate of rja 2 v f 2 ^~L. Comparing the total loss rate to the power input, we find that 



2a 



v f cc 



r)C + a*L/R~ ' (3 ' 6) 
Here £ collects various order-one factors neglected in our rough estimates. We see that 
indeed the front cannot propagate at a constant velocity forever. But v* can stay roughly 
constant until L > (a~ 2 R , and the regime in which this happens is characterized by the 
fact that drag on the membrane behind the front is negligible. (The late-time asymptotics 
have been studied by Granek and Olami |24j].) 

What is the fractional area loss al Suppose first that the modulation saturates behind 
the front at a small amplitude, so that its radius is (Fig. 1) 

r(z,t)=R (l+u{z,t)) (3.7) 

with \u\ <C 1. Requiring that the new shape have the same volume per unit length as 
the old (no global transport of fluid), we find that a must be 0{u 2 ), so that R /a 2 is 
enormous. We will argue that this case is relevant for the case of a short laser pulse, when 
full pearls never form at all. 




Fig. 1: Modulated cylinder. 
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If on the other hand the laser is kept on for a long time, we expect the final state to be 
of large amplitude u ~ 1, so one might worry that a ~ 1. But we know a lot about the final 
state: it must be a string of pearls on thin tethers. After all, neglecting bending stiffness 
the system minimizes area at fixed volume, leading to spherical pearls. Kinetic limitations 
prevent the pearls from all merging; instead we expect a periodic structure with some 
wavelength A. The effect of the curvature stiffness is to prevent the tethers from shrinking 
to zero size. Since a thin cylinder of radius R 1 < R and length I has curvature energy 
^(2irR Q £) (^j-) , Ri will stabilize at a nonzero value. Indeed if the tension is large (as we 
have argued), S/E crit 1 where E crit = k/Rq, then we expect tethers of small radius, as 
seen.i In this case we can neglect the volume and area of the tethers altogether. 

Thus volume conservation says that a string of pearls of radius R 2 , separated by thin 
tethers of length I (so that the wavelength A = 2R 2 + £), must have volume per length 
4?r ^ 2 ' > = 7r(i? ) 2 . The area per unit length is then just A-k(R 2 ) 2 / \. For the wavelengths 
A ~ 2nR characteristic of the observed nonlinear regime, this corresponds to a fractional 
area loss a of just 10%. In other words, the constraint of volume conservation has kept a 
small. If we carry over the estimate leading to ( |3.6| ) (even though now the amplitude is 
not small), we see that once again a~ 2 is very large and since ( is 0(1), we can neglect the 
drag along the length of the tube for the initial propagation. Any gradient of the tension 
must be due to the drag, and moreover we argued in Sect. 3.3 that the tension starts out 
uniform, so the continued smallness of the drag means that the tension stays uniform in 
the initial propagation regime. 

The conclusion of this subsection and the previous one is that in the regime of interest 
we may reasonably model the effect of the laser as generating a constant uniform tension 
on the membrane, just as was assumed in [ r2|]||13|| . Lipid cannot literally disappear locally, 
but in the observed initial regime we may act as though this were the case. In other words, 
the pearling system behaves as if it were suddenly quenched into a uniform unstable state. 
This is the sort of situation in which we expect front propagation, as we recall in Sect. 5 
below. 

As mentioned earlier, our conviction doesn't rest solely on these crude arguments, but 
also on the experimental fact that indeed a front does form and propagate at constant 
velocity, at least initially [12 1 . Later in Sect. 7 we will see the front formation again in the 
numerical solution to a simplified model. 



9 In fact the minimization gives a tether radius of yw/S, times some constants of order 
unity, which is qualitatively correct. There is thus no need to invoke spontaneous curvature 
to explain tether formation. 
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3.5. Other assumptions 

Here we collect other assumptions we are making. 

Our axisymmetric assumption ( j377l) could be generalized to allow u(z, (p, i) = 
Ylim u m( z ' t)e imi ^- The peristaltic mode is m = 0. In Sect. 6.3 we will see that the m > 2 
modes are stable. There remains the m = 1 mode, which corresponds to wandering of the 
cylinder centerline. The m = 1 mode is always soft, or critical, since it corresponds to the 
broken symmetry of displacement normal to the tube. As we mentioned, such wandering 
is indeed observed initially. However we will neglect it since (at least in the linearized 
analysis) it decouples from the interesting peristaltic mode. 

We also assume in our hydrodynamic analysis no-slip boundary conditions between 
the fluid bilayer and the surrounding water. This is a standard assumption and no evidence 
exists to make us relax it. However, the two leaves of the bilayer can slip relative to each 
other, with a friction coefficient measured in the experiments of Evans and collaborators 



[ 21j|||23|l . We will incorporate these effects in Sect. 6 below. 



3.6. Summary 

Let us summarize this long section. We have given a qualitative, intuitive argument 
for a physical picture in which the laser creates a nearly uniform, tense state, both before 
and during the subsequent shape changes. We assumed that the latter proceeded by the 
propagation of a front, which then must initially move at a velocity Vj- considerably smaller 
than the dimensional combination T,/r] (see eqn. ( |3.6| )), as observed. We still need to show 
that this is the case and get a precise formula for Vj-. 

Since front propagation in general leaves behind a state of wavenumber k different 
from the fastest-growing mode /c max [ I9flp0 |) we will also need to revisit our estimate of 



ref. [13] to see what happens to that prediction. We do this for a simple model in Sect. 5, 



then for a more accurate model in Sect. 6. 



4. Very Simple Model 

To get as quickly as possible to the heart of the matter we first develop a simple 
truncated model without a lot of algebra. In this tutorial model only we will make some 
expedient additional assumptions, each of which we will justify or improve in the sequel. 
The assumptions are: 
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i) As in [13], we will temporarily neglect the propagating character of the modula- 
tion until Sect. 5. Instead we will just find the wavenumber of the fastest-growing 
Fourier mode. Our goal here is to show in the simplest possible way how tension 
leads to wavenumber selection. 

ii) We will neglect altogether the bending stiffness of the membrane as well as all 
bilayer effects until Sect. 6.2. This will prove to be a good approximation if the 
applied tension E greatly exceeds the threshold value for instability, as we have 



argued is the case (see eqn. Q3.5D ). 



Hi) We will neglect the finite compressibility of the membrane. This is legitimate as 
long as £ is much smaller than the effective elastic modulus K eff (see eqn. ( |3T5| ) 
and Appendix C). 

iv) We will use lubrication approximation as in (|3.2| ). This is certainly not justified, 
since the wavenumber /c max of the mode we will find is not much smaller than 
unity and hence z derivatives are not much smaller than p derivatives, but it 
will make our equations very simple. Similarly we will neglect the exterior fluid 
altogether. 

We argued that as soon as the laser switches on, very quickly a nearly uniform tension 
£ appears. Volume conservation then implies that the pressure inside jump to a constant 
positive value, p 1 = £/i? (the Laplace pressure of a cylinder under tension) to prevent the 
tube from collapsing. As the tube starts to change shape, described by the small quantity 
u in (|3~7|), the Laplace pressure Ap Laplace (,2, t) = — E • 2H gets new contributions from the 
change in the mean curvature H [|33f : 

^Laplace = ^ [l " « " V Vjti] + 0(u*) • (4-1) 

Since we temporarily neglect the exterior fluid, the interior fluid must move in such a 
way as to give the pressure p{r(z, t), z, t) = Ap Laplace , which we just computed. To linear 
order we again have the fluid flow (|3.2|) , where now /, g are functions of z and time. We 
can no longer fix g by demanding that the net flow of water across any cross-section vanish, 
since now the tube is getting fatter in some places and thinner in others. Similarly we 
need to modify the lipid conservation law. Generalizing ( |3.3|) to variable radius but fixing 
the density gives 

[dz 2Tvr( Zl t)} = -dz-V z [27rr(^, t)v] (4.2) 
15 



where the boundary velocity v z = v e (r(z, £), z, t) = fr 2 — g. Thus linearizing in small 
quantities we get u = —V z [fR 2 — g]. Similarly, conservation of the interior water gives 



d f 

— [dzTvr(z,t) 2 ] = -dz ■ V z (2ivpdp)v 



or«= -V z [\f R o -\9\- Comparing Q we finr* 11 



(4.3) 



u=\R 2 VJ . 



(4.4) 



Since we are linearizing, let us take u(z, t) to be a single Fourier mode: 



u k (t)cos(kz/R ) 



(4.5) 



The fluid equation of motion rj'V 2 v z = V z p, together with ( |4. 1| ) then gives (neglecting 
derivatives compared to p derivatives) 4/ = — ^^-V 2 (l — k 2 )u, or (see Fig. 2) 



u 



8 V R 



k 2 (l-k 2 )u 



(4.6) 



so that the mode with the largest growth rate ii/u is at k max = l/y/2. Note that the tension 
needed to get an appreciable growth rate increases with tubule radius, as observed. 

Note the structure of the dispersion relation ( [4.6|) ; the growth rate is the product of a 
term (oc (1 — k 2 )) arising from the perturbation to the pressure in ( |4.1| ), while the additional 
overall factor of k 2 represents two spatial derivatives coming from the hydrodynamics. The 
first is due to the conservation of fluid volume, and the second from the usual Poiseuille- 
flow relation between fluid flux and pressure gradients. In the more refined hydrodynamic 
theory derived below the general structure is retained, in the sense that the growth rate 
is the product of a hydrodynamic factor and a pressure term, but the former has a more 
complicated wavevector dependence than k 2 as a consequence of the details of flow in 
cylindrical geometry. 



10 Note that in JT3] we mistakenly set the boundary velocity to zero instead of fixing it 
with (^4.2|)-( [13l) . The formula for /c max below is however unaffected. 
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Fig. 2: Exact and lubrication approximation dispersion relations, ne- 
glecting stiffness and bilayer structure. Circles indicate fastest-growing 
modes. 

5. Velocity and Wavelength Selection 

We saw in our very simple model that tension destabilizes a cylinder. More generally 
when we incorporate bending stiffness in Sect. 6 we'll see the same thing, with a finite 
threshold X crit . We now want to know what becomes of such an unstable state. 

We argued in Sect. 3 that in the regime in question our problem could be regarded as 
a uniform, quenched unstable state: the membrane elastically transmits the laser-induced 
tension everywhere. The uniform unstable state changes its conformation until it saturates 
in one of two ways. For brief tweezing, we suppose that the laser pulls out about 1% of the 
area (see Appendix C), then shuts off; even when the laser stops pulling the shape trans- 
formation proceeds to an amplitude such that the projected area is reduced by about 1%, 
i.e. u ~ V0.01. This fits with the observation of a small-amplitude peristaltic shape which 
continues to propagate after shutoff. For sustained tweezing, there is no such limitation, 
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and the shape transformation continues until checked by the curvature energy of the thin 
tethers, as we argued in Sect. 3.4. 

In either case a uniform quenched unstable state rolls off a potential hill to saturate 
nonlinearly. In this class of problems one typically finds that the instability proceeds via 
the propagation of a front from an initial disturbance into the unstable medium. The 
question of the precise conditions under which an equation has such a front solution are 
ticklish and require a more reliable knowledge of the full nonlinear equation than we possess 
(for example, at nonlinear orders the effects of thermal fluctuations are surely more subtle 
than the picture in Appendix C). However, if we take the existence of the front as an 



empirical fact, we can use the analysis of [20] to determine its main properties solely from 
the linearized dynamical equations. 

The shape is governed by some nonlinear, translation-invariant equation in u(z, t) We 
suppose a family of solutions whose envelopes are functions of z — vt, for various values of 
the velocity v. Only one of these solutions can correspond to the actual observed front; 
thus we seek the dynamically selected value v*. The selected value must at least be stable 
to perturbations, in an appropriate sense which we now recall. 

Suppose the front moves from negative to positive z at velocity v* and consider the 
leading edge, where u is still small. We need to examine the effect of small perturbations 
8u. In this region we may use the linearized equations, so that superposition holds and 
the perturbation belongs to a family Su ~ exp(u(q)t — qz) for various values of a complex 
wavenumber q, where oj{q) is the dispersion relation of the linearized equation. The enve- 
lope of this blip is then exp(tReu; — zKeq), which is itself translating at uniform velocity 
v = Reui/Req. Of this family of perturbations, not all are dangerous to the stability of 
the front. A wavenumber q is dangerous only if 

i) Reg > 0, so that this is a leading edge solution, 

ii) Recu > 0, so that it is growing, and 
Hi) v > v*. 

Without (Hi), the blip just gets left behind by the front and never gets a chance to desta- 
bilize it. If any complex q meets all these conditions we will call the front unstable in its 
comoving frame. 



ii 



In the present case the laser itself provides the initial disturbance; for another example 



see the photographs in [38 . 
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Of course the assumed front solution itself corresponds at its leading edge to one 
particular linear solution with some q* , and so 



Reu^/Req* 



(5.1) 



where = u(q*). Let us now consider values of q very close to q* , q 
small a, b. Then we find 



q* + a + ib for 



1 + iw7( R ^-«*)- 



■Imw' 



where u' 



dq 



Since a, b were arbitrary, and q* obeys (i), (ii) above, we can always 



satisfy (Hi) as well unless 



duo 
Im — — 
dq 







„ du> 
Re — 
dg 



(5.2) 



The three equations (|5.1|) - (|5.2j ) are necessary conditions for the linear stability of a prop- 
agating front in its comoving frame |[20|| . We will refer to them as the "marginal stabil- 
ity criterion for v* n or simply "the MSC", and use them to fix the three real variables 
v*, Reg*, Img*. There is a large literature on experimental tests of pattern selection via 
the MSC. Important early studies include refs. |39]. 

We can also use (|5.1|) -( |5.2|) to determine the selected wavenumber q of the stationary 
pattern behind the front. In the comoving frame of the front, the saturated pattern is 
continuously created at the front and moves rigidly with velocity —v*. Nodes are created 
in the leading edge of the front at a rate O = Im(uj^ — q*v*). We can interpret Q as the 
flux of nodes moving toward the saturated pattern; if nodes are not created or destroyed 
as they pass into the non-linear region, then we must have q Q v * = O, or 



q = 1m(u, - q*v*)/v* 



(5.3) 



for the wavenumber of the pattern |1J|. Note that q is generally different from both the 
fastest-growing mode q max and from q* . 

Let us apply (|5.3|) to our very simple model. We found in Sect. 4 that u(k) = 



g^r-/c 2 (l — k 2 ), eqn. (fL6|). Extending this formula to complex values of k and applying 



.20^ ~ 



(pJ] ) - (|5T3|) , we find k = g -^o = '■> a decay length of 3.8i? , and v* 
4 ■ 10 2 / um/sec. In fact we will see in Sect. 6 that our more accurate calculation of k is 
quite close to the very simple model, while our prediction for v* will be much smaller. 
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The selected wavenumber k is purely geometrical since we can form no length scale from 
tension and viscosity; in particular k is independent of laser power. 

Of course the very simple model is still rather crude, even with the refinement we 
just made. We will upgrade it to a more accurate model in Sect. 6. Since the algebra is a 
bit involved, the reader may want to skip this section and pass directly to the simulation 
results in Sect. 7. 



6. Complete Linear Model 

Now that we have some confidence from our very simple model, we need to address 
some of the oversimplifications listed in Sect. 4. First we will replace the lubrication ap- 
proximation with the exact solution to the Navier-Stokes equation. Then we resolve the 
bilayer structure of the membrane to get the four-fluid model mentioned in the Introduc- 
tion. 

6.1. Hydrodynamics: Exact linear theory 

Now we go beyond the lubrication theory and use the full linearized hydrodynamic 
equations for the interior and exterior fluids. We begin with Stokes's observation that the 
axial symmetry we have assumed makes our problem effectively two-dimensional. This, 
along with the incompressibility of the water, 

dv n v n dv 

allows us to write the interior and exterior water velocity in terms of the cylindrical stream 
function ip: 

I dtp I dtp 

v * = -d • (6 - 2) 

The Navier-Stokes equations are now two equations in two unknowns, ip and the pressure 
p. We eliminate the latter by differentiating vjV 2 v p = d p p with respect to z, differentiating 
r]V 2 v z = d z p with respect to p, and subtracting the two resulting equations to obtain 

*-5£)Vo . (6.3) 

In Appendix D we find the stream function i/j required by a single mode of shape 
distortion u(t) expikz/R , enforcing the no-slip boundary condition and lipid conservation. 
This yields v and p in terms of it, and ultimately the stress tensor T+ p inside and outside 
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of the tube which determines the normal force balance equation to relate u and ii, i.e. to 
determine the growth rate u/u = ui{k). As mentioned in the Introduction, our boundary 
conditions differ from |HJ and ||24j| . where the shear stress T zp is taken to be continuous 
across the surface of the cylinder. This boundary condition is appropriate for two viscous 



fluids that meet at an interface |T7[]; in our case, however, there is a material object at the 
interface that supplies whatever forces are necessary to ensure lipid conservation and the 
no-slip conditionii! 

The results for the growth rate have the form that generalizes ( |4.6| ) by the introduction 
of a more complicated "dynamical factor" A(k), 

u{k) = |^A(fc)(l - k 2 ) (6.4) 

with A(k) given by eqn. D.9 of Appendix D. Fig. 3 compares this function with the lubri- 
cation approximation A lub = k 2 /8. Our dynamical factor A(k) differs from the one given 
in |13[ due to the correction mentioned below eqn. (|4.2|) ; however the A(k) of Fig. 3 leads 
to a /c max of .68 (see Fig. 2), very close to the value in [[T3|. 

We see that for all wavenumbers the exact dynamical factor lies below that given 
by lubrication theory. This implies slower growth of instabilities given the same material 
parameters E, k, and rj (see Fig. 2). This increased drag relative to the lubrication ap- 
proximation is sensible; the lubrication model neglects some components of the velocity 
gradients and hence underestimates the rate of viscous dissipation. 

Inserting the growth rate (|6.4[) into the MSC (|5.1| ) -( ^2D and solving the equations 
numerically yields k Q = .75 and v* = .O6E/77. We note that the front velocity is indeed 
much slower than E/77, consistent with our discussion of (|3.6|) . 

We see that taking the propagating character of the instability into account has slightly 
improved the agreement with the experimentally observed k Q . Using our estimated value 
for a, eqn. ( |3.1| ), we also get the front velocity ~ 100 //m/sec. Since our value for E is 



an overestimate we get squarely in the observed range. 



12 In Sect. 6.3 we also address the effect of the two-dimensional membrane viscosity and 
show that it is dominated by the traction of the water. 

21 



0.15 



0.10 



0.05 



0.00 



"i — i — i — i — I — i — i — i — i — I — j — i — i — i — I — r 



/ 



lubrication theory / 

/ 

/ 



0.0 



•3 



complete 
hydrodynamics 



J I I I I I I I I I I I I L 



0.5 



1.0 



1.5 



Fig. 3: Dynamical factor as a function of k 1 with fastest-growing modes 
indicated. 

6.2. Bilayer model 

We now need a more realistic model of the bilayer. For one thing, we have so far 
neglected bending stiffness, which as we saw is eventually crucial to stabilize the tethers 
between pearls. But stiffness effects are comparable to effects involving the differential 
stretching of the two leaves of the bilayer [27] 0, so we must include these as well. This 
means introducing more continuum degrees of freedom, the densities <fi of lipid molecules 
per unit area in each monolayer. 

In fact one must sometimes retain (jr 1 even in equilibrium problems, for example in 
closed systems where the numbers of molecules are separately constrained. This situation 
leads to the "area-difference elasticity" model studied for example in ref. ||. We can 
simply quote their intermediate formula for the elastic energy: 



F[H, <f] =JdS h{2Hf + ^ £ _ i) 2 + E 



(6.5) 
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(For completeness we rederive this formula in Appendix E.) Here k and K Q are the bilayer 
stiffness and compression moduli, and dS is the area element of the midplane. H is the 
mean curvature of the bilayer midplane; in our convention H = —1/2R for a cylinder. 
We measure the lipid density <p + at the neutral surface of the outer monolayer, assumed 
a distance d away from the bilayer midplane, and similarly cp~ . See Fig. 4. Thus the 
parameters of the model are n, K , d, and the full bilayer thickness 2D (see Appendix 
A). Since D and d are much smaller than any other length scale in the problem, we will 



work to leading nontrivial order in them, as indeed we have already done to get ( 6.5|) . 




Fig. 4: Schematic of the bilayer showing various variables. 

We begin in equilibrium at zero tension. Here the densities ^> ± take their preferred 
value (f>Q. Actually we will find it useful to recast our equations in terms of the densities 
referred not to the monolayer neutral surfaces, but to the bilayer midplane. If we imagine 
each molecule casting its shadow on the midplane, the new density variables are the 
excess density fluctuations of those shadows^, relative to the equilibrium densities 0q 
projected to the midplane of the unperturbed surface. Due to the curvature, 0q are no 
longer equal, but instead 4>q = (f> (l^f2H d), where H Q is the equilibrium curvature. Thus 
we have l + x ± =<p ± (lT 2Hd)/<pQ, or 



( j ) ± / ( p = l + X ± ±2(5H)d 



(6.6) 



where 5H = H — H Q is the change of curvature. Substituting in ( |6.5| ) above yields [pi 
F[H,x + ,X~} = J dS{l(2H) 2 + ^[(x + + 2(5H)d) 2 + ( X --2(5H)d) 2 ]+E} , (6.7) 
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In |13[] we called these variables p . 
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Were we to consider fluctuations about flat membranes, we would set H Q = 0; then (|6] 
reduces to the formula in []22[ 



Later we will want to rephrase ( |6.7[ ) in terms of the mean density and density differ- 
ence, and rescale: 



x = (R /d)(x + + x~) , x = (R /d)(x + - x~) 



(6.8) 



We argued in Sect. 3 that the mean density adjusts quickly to its preferred value x = 0, 
since only hydrodynamic drag obstructs this, and initially it is not effective. We cannot 
be quite so sure about density difference, however. The two monolayers are like a pair of 
polymer brushes in contact, and so the friction coefficient between the two layers can be 
(and is) enormous HD 



1. Leaving x f ree 5 we then finally obtain the elastic energy 



F[H,x] = J dS^(2H) 2 + ^f(2H-2H ) 2 



(6.9) 



To gain a physical feeling for (|6.9| ), consider two limiting cases. If the interlayer friction 
is small (or time scales long), we can minimize ( |6.9| ) with respect to x also, and recover 
once again the Canham-Helfrich model. If the friction is very large (or time scales very 
short) then x cannot change at all from its initial value of zero, and we drop the third 
term of (|6.9|) . For flat membranes, H Q = and the first two terms combine to give the 
effective increase of stiffness found in |22 . 



6. 3. Cylinder perturbations 

Let us instead specialize ( |6.9[ ) to the case of a cylindrical initial state and a small 
sinusoidal perturbation of u and x °f wavenumber k. Working to quadratic order in 



fluctuations of shape and density difference one has |33 



H =\R ~ X [-l+u + R 2 V 2 u - u 2 ] 

+ \R [-\{Vu) 2 - 2uV 2 u + (V z u) 2 + 2uV 2 z u] 
dS =(R dzdip)(l + u+ \R 2 {Vu) 2 ) . 
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In these formulas V = (V z , -R 1 ^<p) * s the gradient on the unperturbed cylinder. Substi- 
tuting gives0 

F[u, X] = ^J 2nR dz [const. + (2a - l)u + Q(-R 2 V 2 )u 2 + fx 2 + P{1 ~ k 2 ) X u] , 

(6.10) 

where (3 = K Q d 2 / k and 

Q(x) = l + (a- ±)x + x 2 + (3(1 -x) 2 . (6.11) 

This is our final formula for the elastic energy. As we increase the dimensionless tension 
a (eqn. flOf)) from zero, eventually at some cr crit this quadratic form acquires a negative 
eigenvalue and the system becomes unstable. (In Sect. 4 we dropped all but the a terms.) 

We note in passing that for non-axisymmetric shapes, u(z,(p,t) = u km - cos(kz/R + 
rrup), for m > 1 the polynomial ( |6.11| ) remains stable until larger values of a than for the 
peristaltic m = modes. 

Having found the energy cost F for small distortions in shape and relative density 
difference ( |6.10|) - (|6.11| ), we now need the dynamical equations for the two monolayers, the 
2d fluids of our four-fluid model. These equations are the equations of tangential force 
balance. Each monolayer feels forces due to inhomogeneities in density, the traction T^, of 
the 3d water, 2d viscous forces, and the friction between the two leaves of the bilayer ||22jl : 



SF 

-d zJ - I ±T±+vd 2 d ± Tb(d + -v-)=0 . (6.12) 

Here \x is the 2d lipid viscosity; we have written the frictional force per unit area as 
b(v + — v~) where b is a constant and are the tangential layer velocities ||21|l||23|i [22 



Although the change in x ma Y be neglected just as in Sect. 6.2, now we must keep elastic 
terms corresponding to changes in since these are in principle of the same order as the 
friction term. The second and third terms of ( |6.12| ) are negligible compared to the first 



and the last. Indeed, since v ± ~ u/k, and v + — v ~ -S-ii/k, the last three terms are 



14 This formula differs from the one in |13[] because we no longer need to replace the linear 
terms using the ansatz u = —(u k ) 2 + 2u k cos(kz/R Q ). Instead we will see how the terms 
linear in u determine the constant part of the pressure difference. This rearrangement of 
the algebra has no effect on the answers, but the presentation here emphasizes that global 
volume conservation is not needed and hence the detailed structure of the terminal blobs 
is unimportant. 



25 



of magnitude rju, -^ru, bdii respectively, so that on micron scales the last term dominates 
(see Appendix A). 

To get an equation for d t x, we use the difference in lipid conservation laws d z (v + — 
v~) = —-pj-d t x to write the difference of the dominant terms of (|6.12|) as 



0X + OX' 



or 



k 2 SF 

ot—wi ■ (6 ' 13) 

Next we need the equation of normal force balance analogous to D.7. The Laplace 
force law, eqn. (|4.1|) , needs to be modified to account for bending stiffness. To find the right 
expression we have only to recall the origin of the Laplace law itself: the pressure jump is 
SF/SV, the change in the elastic free energy due to a small change in shape which changes 
the volume by 5V. To get Q we used ff = with ff = R E(1 - R 2 W 2 u) + ■ ■ •) 

and |^ = -Ro 2 (l + u + ■ ■ ■). Now we simply use the full elastic energy ( |6.10| )-( |6.11p . 

We may neglect the difference in monolayer velocities for the purposes of computing 
the hydrodynamic factor A(/c). Once again, this difference is suppressed by d/R relative to 
the central flow velocity and therefore does not significantly modify the energy dissipated 
in the 3d fluids. Thus we again have ii given by —A(k)/r] times the nonconstant part of 
the pressure jump, just as in D.8. Combining with ( |6.13|) , our linearized hydrodynamic 
equations are 



9 fu\ k (A(k) \(Q{k 2 ) + \-a \(3{l-k 



dt\x) Rh\ tk 2 ) { \(3{l-k 2 ) (3 /A 



,2> 



(6.14) 



where e = R rj/bd 2 is a measure of the importance of inter layer friction. The desired 
dispersion relation is then given by the positive eigenvalue, if any, of this linear problem, 
i.e. iv ~ (—180 sec" 1 ) • A, where we took a typical R Q = 0.7/im and A solves A 2 — X(g 1 + 
y) + 9oV = 0- Here g = A(/c) [| — a + k 2 [a — \ ) + k 4 } is the growth rate with zero friction, 
9i = 9o + A(fc)/3(1 — k 2 ) 2 is the growth rate at infinite friction, and y = (3ek 2 /4. 

Fig. 5 shows the growth rate with e = .5, j3 = 3.5, and various values of tension. 
Clearly there is little effect of stiffness on the curve uj{k) when the tension is as large 
as its experimental value a ~ 20 (compare with Fig. 2). For this value we read off a 
fastest growing mode with /c max = .68 for a — 20, but as we argued in Sect. 5 /c max is not 
the expected wavenumber. Instead inserting the above formulas for the dispersion uj{k) 
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Fig. 5: Growth rate versus wavenumber in the full linear model for 
various values of the dimensionless tension a. A typical experimental 
value is a ~ 20. 



into the equations for the MSC ( 5.1 



gives k Q = .80, a decay length of 4.0R , and 



.06^ ~ 



100/mi/sec. The perfect agreement of k Q with Fig. 1 of |T^] is fortuitous, but 
certainly the MSC gives reasonable values for the pattern wavenumber and front velocity, 
and as promised these values are not very different from those found in Sect. 6.1 since at 
large tension both the stiffness and the layer friction are unimportant (see eqn. ( |3.5|) ). 



7. Model for the Nonlinear Regime 

We turn finally to a model for the development of pearls beyond the linear instability. 
A model based on the exact hydrodynamics coupled to the elastic forces presents a rather 
formidable computational challenge. It is therefore of interest to define a simplified model 
which nevertheless maintains certain central features of the full problem. Chief among 
these are (i) conservation of fluid volume within the vesicle, (ii) the monotonic decrease in 
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energy associated with low Reynolds number flow, and (iii) the fully nonlinear structure 
of the elastic energy needed to arrive at a true minimum of the energy functional. 

The simplest model consistent with these constraints utilizes the lubrication approx- 
imation, thus yielding a local evolution equation for the vesicle radius r(z, t). One goal 
of this investigation is to determine whether propagating fronts of constant velocity are 
indeed possible within such a model. Their existence would be strongly suggestive that 
such fronts are supported by the full hydrodynamics. A second goal is to test the validity 
of the marginal stability criterion. While the assumptions of lubrication theory may not be 
quantitatively correct in comparison with experiment, the issue at hand is one of internal 
consistency; within the assumptions of lubrication theory, how accurate is the MSC? 

In the model, the diffusion of tension is ignored, so that E is assumed to be uniform 
along the vesicle. The bilayer bending stiffness is included to stabilize pearls against pinch- 
ing off, but the small effects of membrane compressibility and bilayer friction are omitted. 
Rather than introducing the optical trap explicitly, we consider an initial condition of uni- 
form tension with a localized shape perturbation. We justified this approach at length in 
Sect. 3. 

We first need some formulas for the elastic force beyond linear approximation. The 
pressure difference across the membrane is the change in energy with respect to volume, 
so now our generalization of ( f4.1| ) becomes 

A "-W^r-W- 2 {- EH + K ( t2H + 2H3 - 2HK )} ' <"> 

where V 2 is the covariant laplacian on the curved surface, and K is the Gauss curvature. 
Axially symmetric shapes described by a radius r(z) (Fig. 1) have mean and Gaussian 
curvatures given by 



H = - 
K 



1 / r.. 1 



zz 



2V(l+r2)3/2 r ( 1+r 2)l/2 



r(l + r 2 ) 2 

Likewise, for these shapes the determinant g of the metric tensor and the covariant laplacian 
are 



g = ry/l + r 2 

■ o 1 d r d 



r^l + r 2 dz ^1 + r 2 dz 
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For uniform tubular states of radius R we see from that the balance of forces yields 
an equilibrium radius R = (K/2E) 1 / 2 , as discussed earlier in the context of the thin 
"tethers" between pearls. 

Recalling the derivation of the lubrication theory results in Sect. 4, in particular Eq. 
(fO|) , we know that the equation of motion for the radius r(z, i) should have the form of 
a local conservation law, 

9 2 dJ /~ X 

dl nr =~Tz ' (7 - 2) 

where J is the axial current. For this we appeal to our previous (linearized) result on 
lubrication theory that relates J to the gradient of pressure, 

J=-— r 4 — (7 3) 

4r] dz 

Recall that the factor of 4 in the denominator differs from the usual Poiseuille result due 
to the imposition here of lipid conservation. The pressure is just that determined by the 
energy functional, eqn. (|7.1|) . Combining Eqs. ([7.2]), (|7.3|) , and (|7.1|) , and rescaling both r 
and z with the unperturbed radius R and introducing the rescaled time r = (rjRQ/AKjt 
we obtain the partial differential equation 



dr d 
dr dz 



r 4 — (-aH + V 2 H + 2H 3 - 2HK 

dz V 

where a is again the rescaled tension (|3.1 



(7.4) 



It is easily demonstrated that this is a gradient flow, for the time derivative of the 
functional F is 

dF f 8F dr f , SF 1 / 4 „ 18F 
dz-^— = / dz— cL r o--=— 



dr ./ br dr J 5r r z \ z r 5r 
dir+O. I -— ^1 < 



r 5r J 

where the last line follows by an integration by parts. Thus, provided r > (i.e. the inter- 
face does not pinch off), F is driven strictly downhill. Furthermore, when F is constant 
in time the functional derivative 5F/Sr = and the system is at an energetic extremum. 



Since the model ( |7.4|) also contains the proper linear stability result of the lubrication 
approximation, it provides a dynamics that interpolates between the basic Rayleigh-like 
instability and the stationary final states, while obeying the relevant hydrodynamic con- 
servation laws. 
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The flux form of the equation of motion and the relation between pressure gradients 
and velocity in ( |7.4j) are features found as well in models introduced recently for interface 
motion leading to topology transitions and singularities in viscous flows ||41| . In all of 
these systems, despite the simplicity of the physical ingredients, the equation of motion is 
highly nonlinear and of very high degree (sixth order in ^-derivatives in the present case). 
These features make numerical studies quite delicate, but with care a stable and accurate 
algorithm may be developed. Here we report preliminary results of simulations illustrating 
that propagating fronts of peristalsis are supported by the model. 

Fig. 6 shows the evolution of a cylindrical vesicle with a = 10 perturbed initially 
with a localized distortion. The computational domain has a length 16 x 27r, and periodic 
boundary conditions are imposed by use of a pseudospectral algorithm. The figure clearly 
shows that following an initial "induction" period there is a propagating peristaltic pattern 
moving symmetrically outward from the initial perturbation. A nonlinear least-squares fit 
of the leading edge of the front at successive times to the product of a decaying exponential 
and a cosine function shows that the position of the front increases linearly with time over 
the duration of the simulation, during which it has moved about six pearl lengths. This 
constant velocity is in accord with our physical arguments of Sects. 3-4. We find the 
dimensionless wavenumber of the pattern in the leading edge of the front to be k ~ 0.78, 
and the decay length to be ~ 5.4. Applying the techniques of the earlier sections to the 
model of this section, one finds that the MSC predicts a wavenumber k Q = 0.77; the decay 
length of 4.3 is somewhat smaller than the value observed in the numerical solution, but is 
still on the order of a pearl size. Thus the front is "sharp." Finally, in the rescaled units of 
([7.4|) , the front velocity at the value a = 10 is found to be ~ 13, close to but somewhat 
smaller than the marginal stability criterion value of 17. Work is in progress to understand 
the origins of these discrepancies. 

In addition to these quantitative aspects of front propagation, we see that the state 
left behind the front indeed has the appearance of a string of pearls, with very narrow 
tethers between somewhat prolate ellipsoids. These are quite similar to ones seen in the 
experiments of Bar-Ziv and Moses. It is also interesting to note that a slight shift between 
the selected wavelength near the leading edge of the front and that found well behind, 
where the pearls are fully formed. This difference is in qualitative accord with what one 



finds in other more familiar examples of propagating pattern selection |L9 
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Fig. 6: Numerical solution of the lubrication equation (|7.4j ) with a = 10, 
showing the formation of pearls. Time increases from bottom to top in 
increments of At = 1.0. 

8. Conclusion 

Laser tweezers are an important new experimental tool for understanding membrane 
dynamics; they furnish a precisely controllable, local, physical intervention. In this paper 
and |y§ we have proposed a simple model for the effective laser-membrane interaction and 
used it to explain much of the phenomenology of the pearling instability discovered in this 
context by Bar-Ziv and Moses. 

Various aspects of the pearling phenomenon remain to be explained. For example, 
we have not attempted to study the migration of the pearls which develop after prolonged 
tweezing. What we did show was how a cylindrical vesicle under tension is unstable to 
states with a periodic modulation in diameter. We refined our previous prediction for 
the initial wavelength of this perturbation using the marginal stability criterion and pre- 
dicted a selected wavenumber in good agreement with the experiments of ref. [|T2|. We 
also explained the observed front behavior and predicted a front velocity in qualitative 
agreement with the experiments. Along the way we showed how to incorporate the bilayer 
membrane structure into the dynamics. Although this structure did not significantly affect 
our results, it may matter in dynamics problems at lower tension and on shorter length 
scales. Finally, our numerical calculations have demonstrated that the lubrication approx- 
imation is rich enough to capture the essential physics of the pearling phenomena, such 
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as the existence of an initial uniform front propagation velocity and the shape of the fully 
developed pearls. These calculations were in reasonable agreement with the predictions of 
the MSC, thus bolstering our faith in this approach and suggesting that the MSC can be 
reasonably applied to other fluid dynamics problems. One particularly relevant problem 
is the propagation of the ordinary Rayleigh instability along a column of fluid, as seen in 
recent experiments on droplet fission [[12]. 

We are not aware of any naturally-occurring biological process to which the pearling 
instability literally applies. Our goal was rather to create a simple theory of membrane 
dynamics and test its many underlying assumptions by applying it to a well-controlled 
dynamical shape transformation. We expect that the approach developed in this paper 
can provide a starting point for understanding other phenomena of more direct biological 
interest, for instance the dynamics of shape instabilities during adhesion and budding. 
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Appendix A. A Compendium of Constants 

Here we list the values of various physical constants for DMPC bilayers. We take: 
the bilayer bending stiffness to be k = 0.6 • 10 _12 erg 



the bilayer thickness to be 2D = 40 A [44 
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the distance from bilayer midplane to each monolayer neutral surface to be d ~ D/3 = 12 A 

Hi, 



the bare bilayer compression modulus to be K = 144 erg cm -2 f4"3 
the interlayer friction constant b to be 10 8 ergseccm -4 pTf , 
the viscosity of water r\ to be 1 • 10 _2 ergseccm -3 , 



and the 2d viscosity of the lipid layer to be fx = 10~ 6 erg sec cm" 2 [f4"5f . 

In addition for concreteness we sometimes take: 

the typical tubule radius R = 0.7/xm(ref. |i~2|| , fig 1), 

the typical laser power 50mW fi~2Tl , 

and the temperature to be k B T = 4.2 • 10 _14 erg . 

Thus we get the derived quantities: 

k B T/4nK = 0.006, 

E crit = k/R 2 = 1.2 ■ 10- 4 ergcm- 2 , 
e = R Q r]/bd 2 ~ 0.5, 
[3 = K Q d 2 /K~3.5. 

Also, we estimated the induced tension E ~ 2 • 10 _3 ergcm -2 in Sect. 3.2, so 
Yjjr\ ~ 2 • 10 3 /xm/sec and the dimensionless tension 
a = ZR 2 /k ~ 20. 

The values of E, b, d are not known very accurately, but fortunately our formula for k is 
not very sensitive to them. 



Appendix B. The Delaunay Surfaces 



We argued in the text against the introduction of a spontaneous curvature term in 
the elastic free energy. In the presence of such a term we may add and subtract a constant 
to recast the energy in the form 

F[shape] =2k [ \(H - H ) 2 + a' 



Large thermal fluctuations then imply that initially a' = 0, and the free energy then simply 
prefers conformations with constant mean curvature everywhere equal to H . 

Remarkably there is a one-parameter family of such shapes with axial symmetry. 
These are the "unduloid" or "Delaunay" surfaces. They obey a nonlinear equation ex- 



pressing H = H Q . Deuling and Helfrich used this equation in ref. ||32|| . There is however a 



very simple geometrical characterization of these surfaces [46 
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Draw a line I in the plane. Construct any ellipse, tangent at one point to £. Let 
the ellipse roll without slipping along t, and follow the path of one focus to get 
a curve C. The Ggure of revolution obtained by rotating C about the axis I is a 
surface of constant mean curvature l/2r , , where r + is the larger semimajor axis 
of the ellipse. 

Let us see what this says for nearly-cylindrical shapes. Such shapes are generated 
from nearly-circular ellipses, with foci near the center. The radius R of the near-cylinder 
is thus equal to that of the near-circle. Rolling the near-circle along the axis of course 
gives a periodic shape with period A = 2ttR q . In our language, this is a small perturbation 
with wavenumber k = 1, larger than the initial wavenumber seen experimentally. 



Appendix C. Thermal Fluctuations 

In subsection 3.3 we quoted an estimate for the effective stretching modulus of the 
bilayer. 

Generally the effects of thermal fluctuations on the conformations of stiff membranes 
are quite small unless one carefully adjusts to the threshold of a shape transformation. One 
simply replaces the elastic constants by effective values corrected by terms of relative order 
* B T/4™, about half a percent for DMPC.B We wiH neglect such effect, One notable 
exception, however, is the susceptibility K , whose "bare" value is so small that thermal 
corrections can be significant. At modest values of tension, the membrane effectively 
becomes elastic, because significant area can hide in invisible short-scale wrinkles excited 
by thermal agitation [[7] . (The tiny value reported for the bare Kq 1 in is attained only 



at extremely large values of applied tension, where these wrinkles have been extinguished.) 

Following Helfrich and Servuss, let us estimate the magnitude of this effect in our case. 
We consider a membrane thermally fluctuating about a cylindrical shape, and evaluate the 
average true area A = (f dS) in terms of the "projected" or apparent area 2nR Q L. We 
would like to know how A changes as we adjust the tension, at fixed projected area; this 
tells us how much area we can pull out of our membrane at fixed apparent shape, and 



hence the effective modulus. Following | 48|j , the relative area extension as we change the 
dimensionless tension a is 

.... k B T f^dkk 2 



8tt 2 J.^ 



15 Something similar happens with dynamic couplings 
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In our case a ~ 20 (Appendix A.). The lower bound of the integral should be either the 
initial tension or unity, whichever is larger, to account for the cylindrical geometry; since 
we take the initial tension to vanish and a ~ 20 we get 6 A/ A ~ (k B T/8nK) log 20 ~1%. 
This estimate is not very sensitive to the actual value of a. Then the effective modulus is 
given by the constitutive relation E = K eS (5A/A), or K eS ~ 100 • E ~ 2 • 10 _1 ergcm -2 , 
so that the hierarchy (|3]5|) we assumed earlier is still preserved. 



Appendix D. Computation of the Dynamical Factor A(k) 

Eq. ( fOp is solved by separating variables as if) = ^(p) exp(ut + ikz/R ). ^ satisfies 
(dp — p~ l d p — (k / Rq) 2 ) 2 ^ = 0, which has solutions of the form 

= A.pl, (kp/R ) + B x pK x (kp/R ) + A 2 p\ (kp/R ) + B 2 p 2 K (kp/R ) , (D.l) 



where I u (x),K u (x) are modified Bessel functions | 49fl . Demanding that ^(p) be well- 
behaved at p = and p = oo forces B l = B 2 = in the interior fluid and A x = A 2 = in 
the exterior fluid. The no-slip boundary conditions, together with the lipid conservation 



law (|4.2|) , give us four equations for the four unknowns A x , A 2 , B x , B 2 : 

v+(p = R )=R u, v p -(p = R )=R u, (D.2) 

vt(p = R )=v = ^u , v -{p = R Q )=v= T ^u , (D.3) 

where + (— ) denotes the exterior (interior) fluid. Note that in ( p.3| ) we have assumed 
that the lipid velocity v has no constant part, in accord with the arguments of Sect. 3 that 
after the tension has spread the 2d pressure (tension) is practically uniform and we can 
take any net flow of lipid to be small. Also, we have taken the lipid to be incompressible 
since the applied tension is much smaller than the membrane compressibility (eqn. (|3.5|) ). 
Plugging ( p.l|) into ( p.2|) - (|D~3|) , and using the relations 



I' (x) = l 1 (x) , l' 1 (x)=l (x)-h 1 (x) , (D.4) 

K' (x) = -Kj(x) , K[(x) = -K (x) - iKj(x) , (D.5) 
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we find 

i^-( A i l i(k) + A 2 R l (k)) =R ii 

i^-{B 1 K 1 {k)+B 2 R K (k)) =R Q ii 

A^-W + A 2 (2I (fc) + k\(k)) =^u 

B 1 ^-K (k) + B 2 (kK 1 (k)-2K (k))=^u , (D.6) 

i.e. the A's and the i?'s are proportional to ii, with proportionality constants which are 
rational functions of the modified Bessel functions. 

To relate ii to u and thus determine the growth rate u/u, we examine the normal force 
balance equation at the membrane: 

Tp P ~ T~ - Ap Laplace = , (D.7) 

where is the 3d stress tensor of water and — Ap Laplace is the pressure jump due to the 
surface tension, given by (^4.1| ). The z component of the Navier-Stokes equation relates the 
pressure to the velocity, p ± = ^^-V 2 ^ + Pq ■ The NS equation does not fix the constants 
Pq] we must choose p, = and p, = £/-R . Using 

for f(p) = I (^), K (Jg) ||, we find at the membrane^ 

T+ ~ T- = - 2ir,^-(B 1 ^-(K (k) + l Kl (*)) + B 2 kK,(k)) 

- X^A^Qvih) - ^(k)) + A 2 k\{k)) + E/R (D.8a) 
= -r]A(ky 1 u + 'E/R . (D.86) 

Eqn. ( P-8EQ defines what we call the "dynamical factor" A(/c) in Eq. ( |6.4|) . Taking the 
values for the A's and S's determined by ( |D.6| ), we find 



Art x 1 [k(Kj - Kf) + 2K K X ] [k(I$ - If) - 21,1,] 
2k 3I K /k - KJJk + k{lfKl - llKf^ 



In eqn. ([D.9|) all Bessel functions are to be evaluated at k. 



16 We have used = —pS^ + ^(c^u ■ + dj v i) f° r incompressible fluids |5(| 
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Appendix E. Membrane Elasticity 

Here we briefly derive eqn. (|6.5|) . Following previous authors (e.g. HU) we visualize 
each monolayer as a fluid of compressible cylinders, each of length D. One end (the "chain" ) 
touches the bilayer midplane and resists compression and extension with a modulus K c 
and preferred area/molecule at the end of a c0 . The other end (the "head") is a normal 
distance ±D from the midplane with modulus K h and preferred area a h0 . All of the 
parameters introduced are the same for each monolayer, since the two layers and their 
respective solvents are identical. We thus consider a monolayer elastic energy per molecule 
of 

/ ± \ 2 / ± \ 2 

/ ± [shape,aJ,a c ± ] = i^a J^-l + \K c a 0c ^ - 1 -const. . (E.l) 

\ a 0h J \ a 0c / 

Here the constant includes the chemical potential for extracting additional molecules from 
any reservoir; it will of course turn out to be related to the physical tension. 

It's convenient to introduce a pair of imaginary "neutral" surfaces at normal distance 
±d from the midplane, where d will be specified in a moment, and to measure the densities 
(f)^ of lipid molecules in the two monolayers relative to these surfaces. Thus the areas c 
above are 

4 = (<P ± r 1 [lT2H±(D-d)} 
af = (</) ± )~ 1 [l ±2H ± d] . 

Here #± are the mean curvatures of the two neutral surfaces. We dropped terms involving 
Gauss curvature, as well as those of order (HD) 2 . Substituting in ( |E1| ) we get 



const.' 



3 



(E.2) 



where K, <fi , k, c m , the new constant, and the neutral surface location d are various com- 
binations of the original parameters chosen to put ( |lil.l| ) into the form ( |E2| ) . In particular 
d is chosen to set the coefficient of the cross-term in the convenient form shown. Since the 
original parameters weren't directly observable anyway, we abandon them now in favor of 
the new ones. This is actually progress. Note that again the new parameters are the same 
for each monolayer. 
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The monolayer free energy is now J dS± where dS^ are the areas of the two 

neutral surfaces. Rewriting in terms of the area element dS of the midplane we get 



dS ± (f^r = dS(lT2Hd) 



#0 



(^^^C^ + ^T -const/' 



Here if is the mean curvature of the midplane: if ± = H ± 2H 2 d plus a Gauss-curvature 
term we drop. 

We will be considering only values of the tension much smaller than K (eqn. ( |3.5| )), 
so we may use </> in place of 4>^ in the constant term above. Writing the other H in 
terms of H and finally identifying the constant with the tension E then yields eqn. (|6.5|) . 

Before switching on the laser the terminal blobs serve as reservoirs, and we argued in 
the text that then E should be taken to be close to zero. In equilibrium the densities will 
then just take their preferred values, and we recover the Canham-Helfrich model, with k 
the bilayer stiffness. In particular, no spontaneous curvature appears, as argued in pi[| . 



After switching on the laser the system is no longer in equilibrium, and we need to work 
out dynamically what happens. 

We note in passing that in a closed system the constant term in ( |E.1|) will not be equal 
for the two layers; rather than containing a chemical potential for a reservoir, these terms 
act as La grange multipliers enforcing the constraint of fixed molecule numbers in each 
layer. Letting = a Q N^ be the "preferred area" of each layer, we then recover the 
area-difference elasticity model (ref. 0). In particular the parameter ct ADE is essentially 
the same as our parameter f3. 
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